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An analytic approach to phenomenological models inspired by cubic string field theory is in- 
troduced and applied to some examples. We study a class of actions for a minimally coupled, 
homogeneous scalar field whose energy density contains infinitely many time derivatives. These 
nonlocal systems are systematically localized and an algorithm to find cosmological solutions of the 
dynamical equations is provided. Our formalism is able to define the nonlocal field in regions of the 
parameter space which are inaccessible by standard methods. Also, problems related to nonlocality 
are reinterpreted under a novel perspective and naturally overcome. We consider phenomenological 
models living on a Friedmann-Robertson-Walker background with power-law scale factor, both in 
four dimensions and on a high-energy braneworld. The quest for solutions unravels general features 
of nonlocal dynamics indicating several future directions of investigation. 
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I. INTRODUCTION 

The interest in nonperturbative formulations of string 
theory has been increasing for two complementary rea- 
sons. From a theoretical point of view, the dynamics of 
strongly interacting strings and its consequences in the 
articulation of a fundamental, predictive theory are far 
from being fully explored. In parallel, cosmology is a 
natural field for testing string theory, since observations 
at very large scales (e.g., of the cosmic microwave back- 
ground) probe energies close to the Planck scale. As soon 
as new developments are achieved on the theoretical side, 
phenomenological details are gathered on their cosmolog- 
ical aspects. Ultimately, we hope that this process is two 
way and cosmology will be able to indicate promising 
research directions to the string community. 

One framework falling into the above discussion is 
string field theory, particularly in its cubic version 
(CSFT); see [1 S (3 for some reviews. The effective 
action of the open CSFT particle modes is nonlocal, 
in the sense that the fields are derived infinitely many 
times by a set of d'Alembertian operators of the form 
e r»n _ 2 _|_ r j2 + ?'^D 2 /2 + . . . , where r* is a number 
fixed by the theory. The simplest (and best studied) ex- 
ample of such effective actions is obtained when all fields 
are switched off except for the scalar tachyon. In the 



Minkowski case, one can find that approximate (although 
very accurate) solutions satisfy a diffusion equation with 
respect to an auxiliary evolution variable r, eventually 
set equal to the physical value r* 0, [j| . 

One can try to apply this procedure also in the pres- 
ence of gravity. There are a number of reasons why 
to consider the cubic SFT tachyon also on a cosmologi- 
cal [Friedmann-Robertson-Walker (FRW)] background. 
First, nonlocal theories are still poorly understood and 
it would be instructive to embed them in an expanding 
universe and compare the resulting cosmology with the 
standard big bang/inflationary model. Second, both the 
bosonic and supersymmetric Minkowski solutions found 
in [1 [H have a spike at the origin of time. 1 One might 
ask whether a Hubble friction term would modify this 
behaviour in the FRW extension of such solutions. 

Third, by now a cosmological tachyon has been exten- 
sively studied in its Dirac-Born-Infeld (DBI) formula- 
tion. The first version of the DBI effective action [6| (see 
0, H| and references therein for an overview of the liter- 
ature) is valid when higher-than-first derivatives of the 
field are small compared to its speed d, EE EH, EH ■ The 
simplest versions of tachyon inflation do not reheat and 
generate large density perturbations, unless the param- 
eters of the model are fine-tuned [11 El El EE 53, EE 

tHE I1H- Although these problems can be bypassed 
, it is unlikely that even a viable tachyon would lead 
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1 Unless one promotes the regulator in the solutions as a physical 
parameter; in this case, the scalar field is a regular function of 
time. 
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to predictions qualitatively different from those of stan- 
dard inflation [8j. Moreover, DBI actions tend to de- 
velop caustics on an inhomogeneous background, which 
break down the DBI approximation (higher derivatives 
diverge on caustics) [23|]. It was also shown that the 
field is unsatisfactory as model of quintessence [HI - 2 In 
a more recent and promising formulation of DBI scenar- 
ios, higher-order derivatives are taken into account and 
new features, such as large non-Gaussianities^can emerge 

m m m m, m m m m m m m m Mm. 

There is evidence that the cubic SFT tachyon is not 
plagued by the above fine-tu ning s and m ay g ive interest- 



ing phenomenology [aa [40|, [4JJ, [42, 1431 [4^, [451, [46|, J47 

IH,|4§|. However, there are several conceptual issues 41 1 
which have not yet been addressed. It is our purpose to 
verify these issues and comment on past results. We start 
this research program from the study of phenomcnolog- 
ical models, presenting the detailed CSFT tachyon case 
in a forthcoming publication (Hoj . 

It will be convenient to inspect the effective cosmol- 
ogy on Randall-Sundrum (RS) and Gauss-Bonnet (GB) 
four-dimensional braneworlds, in the approximation of 
hi gh energy and no brane-bulk energy exchange (see 
[& l4ll] for an overview of the braneworld formalism used 
here 3 ). The local equations are easily generalized to these 
cases [H, which are interesting in their own because the 
superstring tachyon mode naturally lives on a non-BPS 
brane @. 

It is important to stress that the methodology of this 
paper and its companions 0, which has tight con- 
nections with the 1 + 1 Hamiltonian formalism developed 
in [H, [H, Hi], [H, [56| , is of particular relevance when the 
scalar field is interpreted as the inflaton. There, the field 
is quantized and perturbed to compute its n-point corre- 
lation functions (primordial (n — l)-spectra), an aspect 
which will be not dealt with here. Without localizing 
the derivative operators, it is difficult to define a set of 
fields, conjugate momenta and commutation relations in 
the usual Ostrogradski way. To the best of our knowl- 
edge, there are no other known ways to achieve a sensible 
quantization while retaining the genuine nonlocal nature 
of the problem. 

Another benefit of our proposal is that the nonlocal 
Cauchy problem [57| becomes transparent and reduces 
to the determination of a finite set of parameters of the 
candidate solution. The reader is referred to Eqs. (|86|) 
for a braneworld example. 
This is the outline of the paper: Section |TT] contains 



2 In general, the string tachyon is regarded as a late-time candidate 
for dark energy only on phenomenological grounds, because its 
evolution scale is of order of the string scale. Tachyon dynamics 
can be associated more naturally to the very early universe. 

3 This approximation at the level of the background equations was 
invoked also to simplify the analysis of large-scale inflationary 
perturbations. The legitimacy of this attitude has been recently 
validated by the results of the full computation of Ref. [5lll . 



some review material on cubic SFT and nonlocality on 
curved (in particular, FRW) backgrounds. In Sec. |TTT] ex- 
act solutions of the local system are reviewed. Section Hvl 
describes an attempt to find nontrivial cosmological so- 
lutions via a covariant Ansatz, while in Sec.|V]it is shown 
that nonlocal solutions as infinite power series cannot be 
constructed from exact local cosmological solutions. Sec- 
tion |VI] explains the strategy to find approximated nonlo- 
cal solutions with the local solutions as initial conditions 
in an auxiliary evolution variable. These steps are ap- 
plied concretely in Sec. I VIII for braneworld cosmologies 
with constant slow-roll parameters; there, we confront 
the truncated power series of Sec.[V]with the "localized" 
function thus found. The same procedure is applied in 
Sec. I Villi for 4D power-law cosmology. A discussion of 
the results, comparison with literature, and open prob- 
lems are in Sec. IIXI The Appendix contains technical 
material. 



II. SETUP 



The bosonic cubic SFT action is of Chern-Simons type 



S= - 

9, 



(1) 



where g is the open string coupling constant (with 
[g 2 ] = E 6 ~ D in D dimensions), J is the path integral 
over matter and ghost fields, Qb is the Becchi-Rouet- 
Stora-Tyutin operator, * is a noncommutative product, 
and the string field $ is a linear superposition of states 
whose coefficients correspond to the particle fields of the 
string spectrum. 

At the lowest truncation level [59(, all the excitations 
in <3> are neglected except for the tachyon field, labelled 
4>{x) and depending on the center-of-mass coordinate x 
of the string. The Fock-space expansion of the string 
field is truncated so that <E> = |$) = <j>(x)\ J.), where the 
first step indicates the state- vertex operator isomorphism 
and I J.) is the ghost vacuum with ghost number —1/2. 
At level (0, 0) the action becomes, in D = 26 dimensions 
and with metric signature ( — h . .+) [60L l6lj|. 



S4, - -T 
So 



(V 



^(a'd^ + 1)0 

3 



A 



(2) 



where A = 3 9 / 2 /2 6 w 2.19, a' is the Regge slope, and 
Greek indices run from to D — 1 and are raised and 
lowered via the Minkowski metric rj^ . The tachyon field 
is a real scalar with dimension [(f)] — E 2 . The constant 
A does not contribute to the scalar equation of motion 
but it does determine the energy level of the field. In 
particular, it corresponds to the D-brane tension which 
sets the height of the tachyon potential at the (closed- 
string vacuum) minimum to zero. This happens when 



3 



A = (6A 2 ) -1 , which is around 68% of the brane tension; 
this value is lifted up when taking into account higher- 
level fields in the truncation scheme. 

On a curved background, we define the operator (in a 1 
units) 



A°/ 3 = 



.□_vSG*A)' 



E 

1=0 



□< = $>□< 

1=0 



where 



and 



□ 



r» = — = ex = ln3 3/2 - In 4 w 0.2616. 
3 

Defining the "dressed" scalar field 

and coupling the tachyon to a general metric 
total action is 



(3) 



(4) 



(5) 



(6) 



the 



S 



R 

2k 2 



1 



+ - </>(□ - m )<f> — U(<j>) — A 



where g is the determinant of the metric, R is the space- 
time Ricci scalar, k is the effective gravitational coupling, 
to 2 is the squared mass of the scalar field (negative for the 
tachyon), U is a nonlocal term that generalizes the cu- 
bic potential of Eq. ([2]) , and we have absorbed the open 
string coupling into (j). For simplicity we have written 
the gravitational part as the Einstein-Hubert action, al- 
though later on we shall discuss other possibilities. 

We specialize to a four-dimensional target where the 
other dimensions are assumed to be compactified, sta- 
bilized, and integrated out. On the flat FRW metric 
ds 2 = —dt 2 + a 2 (t) dxidx 1 , for a spatially homogeneous 
field (j> = (j>(t) the d'Alembertian operator reduces to 



□ = -(d 2 + 3Hd t ) . 



(8) 



where t is synchronous time, a(t) is the scale factor, H = 
a/ a is the Hubble parameter, and dt = ' is the derivative 
with respect to t. All the formulae below reduce to the 
Minkowski case for H = 0. 

Following the notation of [41] (to which we refer for 
the detailed derivation of the dynamical equations), the 
equation of motion for the nonlocal scalar field is 



D4> = m 2 <j) + U' . 



where 



U 1 



U' 



dU 



(9) 



(10) 



is constructed from a nonlocal potential term U (4>) which 
does not contain derivatives of <p. One can also recast 
Eq. ([9]) in terms of </>, 



(□ - m 2 )e 



= U'. 



(11) 



In the string cases the potential is monomial, U(<f>) = 
<j) n , where er(r*) is a coupling constant (n = 3 for 
the bosonic open string, n — 4 for the approximate su- 
persymmetric string). The total potential is 



v(<M) = -TO 2 <t? 



u(4>) + A . 



(12) 



On a FRW background there is an extra equation relating 
the gravitational sector to the scalar one: 



/frw(^) 

where 

0i 



H 2 ~ — 
3 



Y^ 1 



o 2 ) + v-o 1 



= 0, (13) 



r* / ds(e sr ' a U')(ne 



r* 



f ds[{U-m 2 )e {s ~ 2)r * U 4) 
Jo 



<P](ne- sr * a $), (14) 



02 = 



2r» 



dsd t {e sr * a U')d t {e- sr * a (, 



dsd t [(D 



m 2 )e ( s -2)r.D 



4>}d t (e 



-sr# □ ' 



(15) 



In the local case (r* = 0, A = 1), Oi = 0. 

It is straightforward to generalize Eq. (TIB"]) to scenarios 
where the visible universe lives on a (3+l)-brane embed- 
ded in a higher-dimensional spacetime. If one neglects 
energy exchange between brane and bulk modes, the only 
modification to the equations of motion is in the Fried- 
mann equation (|13[) . which becomes, at high energies, 



0. 



(16) 



The constants 9 and /3g are determined by the brane 
setup. In particular, for a Randall-Sundrum brane 9 = 1 
[H, IH, [6J, [65|, J6f| . while for a brane in Gauss-Bonnet 
gravity 9 = — 1 [67], [H, [6§| . At low energies, one recovers 
general relativity. We shall comment on the approxima- 
tions implied by Eq. (|16p in the final section. 



III. EXACT LOCAL SOLUTIONS 

Before considering the nonlocal system, we recall some 
exact solutions of the standard local equations of motion 
(r* = 0, <j) = 0) with to 2 = (mass terms can reappear 
in U). We start from a configuration 



H 



<pot p , 
H t q , 



(17) 
(18) 
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where p and q are real numbers and (f>o, Hq are constants. 
This Ansatz parametrizes two widely studied classes of 
exact cosmological solutions of the local equations, both 
in general relativity and other theories of gravity such as 
the braneworld. 

The first class has a power-law scale factor (q = — 1, 
a = t Ha ), corresponding, in four dimensions (9 = 0) and 
modulo integration constants, to [70j, l7l| 



2H 



In t . 



(3H - 1)H 
U = ~ exp 



(19) 
(20) 



In the Randall-Sundrum (6=1) and Gauss-Bonnet 
— 1) braneworlds, 



kt o/2 = lJlzl H a-o)/^e/2 



e V 3/2 



(21) 



u = li 



2-6 / ± \ 2(1-2/0) 



QH J 0? 



e \<P0 



A 2(l-2/0) 



(22) 



which are 4> t 1 ^ 2 1 U oc 4>~ 2 and <\> oc i^ 1 / 2 , U oc 6 , 
respectively (see [72( for details). 
Another exact solution reads 

a(t) = exp(pt n ) , H = pnf 1 ^ 1 , sgn(p) = sgn(ri) , 

(23) 



which gives 



= o t [n(1 " e)+e]/2 , 

2(ti-1)(2-9)-2ti 

U = A<f> «(!-«)+« 



2(rt-l)(2-9) 
i n(l-fl)+e 



(24) 
(25) 



for some A = ^(n,p) and £ = B(n,p) [72j,l73j,l7|. There 

is also a Gauss-Bonnet solution cx In t with exponential 

potential when n = 1/2 [72], which we shall not consider 

here. Finally, the case n = 2 (linear Hubble parameter) 

may be of particular interest because in four dimensions 

it corresponds to a quadratic potential: = (jyot 1 ^ 8 ^ 2 , 
U = A(j) 2e/(9-2) +B ^2_ 

We note that it is possible to formally get a potential of 
type Acj) 2 + B(j) 4 , A ^ ^ B, only in two situations. The 
first is for 9 = in the limit n — * oo in Eq. (|25[) . where 
the coefficients of the potential as well as a and -ff are 
ill-defined. The second is an effective scenario (with no 
connection to any known modified gravity model) with 
9 = 4 and n = 2, which has <f> = <po/t and a Hubble 
parameter linear in t. The effective potential is then 



IV. NONLOCAL SOLUTIONS: COVARIANT 
ANSATZ 

The simplest way to solve the equations of motion is 
by requiring that 



□</> = U + A<j), 



(27) 



where Uq and A are constants. This Ansatz, first intro- 
duced in [75[ for nonlocal gravitational degrees of free- 
dom, has the advantage of being covariant, and hence 
independent from any specific background. It does not 
require knowing the Hubble parameter H (t) to compute 
any nonlocal function like the CVs. The aim of this ex- 
ercise is twofold. On one hand, it will show that there 
exist solutions of the nonlocal equations which are trivial, 
inasmuch as they behave just like local solutions of local 
models with rescaled parameters. On the other hand, by 
following Eq. (|2"?| we will be forced to abandon string 
theory and limit ourselves to a phenomenological model, 
letting the field mass be a free parameter and the nonlo- 
cal potential be unspecified. The simplifications given by 
Eq. P?| are remarkable, 4 but globally incompatible with 
string models. The key point is that the eigenvalue equa- 
tion of the d'Alembertian operator is still the right path 
towards an analytic treatment of nonlocal systems, but 
under a different interpretation than that of Eq. (12~T1) . As 
we shall see, the eigenfunction of the □ will not be the so- 
lution itself, but the kernel of an integral representation 
of the latter. 

Comparing Eqs. ([9]) and (|2T[) one has that A = m 2 and 
U' = U a . Also, noticing that U n <j) = (U + m 2 0)(m 2 )"~ 1 , 
one gets 



U (l 
U + 



) + m e 



2r,m^ 



Ucj>, 



U4> = m 2 {4>-U ) 



(28) 
(29) 



The last equation and Eq. (TTTj) imply that U' = U' = 
B(Uo,m 2 ), U = B(Uo,m 2 )(f> + C, where C is an integra- 
tion constant and 

B(U ,m 2 ) = ~m 2 U (e~ 2r * m2 + 1 + e r * m2 ). (30) 

Hence, 

Ox = B(Uo,m 2 )(e r * m ' 2 - 1)(U + m 2 0), 2 = 0, 

(31) 



U = ip{4> 2 - 2(5\p 2 ^) 



(26) 



which has an extra — sign relative to the string potential 
if the universe expands and 4> is real (p > 0, /3| > 0). 
This has the shape of a (p = 3)-adic potential. 



4 With a quadratic potential, the rescaling of the parameters can 
be used to heal the r\ problem and obtain scale-invariant pertur- 
bation spectra also from local potentials which are usually too 
steep to support slow rolling. The effective parameters of the 
model are chosen to tune the slope of the potential while the 
physical ones maintain natural values [4§|. 
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so that the system reduces to a standard, local chaotic 
inflationary model with energy density 



D(Ub,m 2 ) = B(^ ,m 2 )[l + (l 



-)m 2 ], 



A(C,C/ ,m 2 ) = C + S(i7 ,m 2 )(l-e r - m )U . (32) 

One can eventually choose C so that A = 0. This system 
has three free parameters {Uo,m 2 ,C} as the standard 
phenomenological model. 

An example of a trivially local asymptotic solution is 
de Sitter evolution with H = s/C and <j> cx e _t . More 
generally, one can assume an exponential Ansatz in de 
Sitter (if = H ): 



= 2J a n ( 



b n i 



(33) 



A. de Sitter case (q — 0), 4> = t p 



Although an exactly constant Hubble parameter re- 
quires the scalar field to be nondynamical, the dS solu- 
tion is a good description of the standard inflationary 
regime. 

In order to find <j>, one has to act with multiple □ 
operators on the scalar field. It can be shown that 



.t-n 



n+e-i 

n (p-v 

T{£ - p)V(-£, 1-21 + p; -3H t) 



n=0 



n-p) 



±p—n—i 

i-iytp- 2i , 

(37) 



so that 



E 



-r,b n (b n +3H )t 



and the scalar equation of motion reads 
- 2J a n [b n (b n + 3H ) 



m 2 ]e Kt 



(34) 

where $ is the confluent hypergeometric function of the 
second kind described in the Appendix. For fixed t and 
Hq, ^ — > 1 when £ — > 0. The series defining </> is not 
pointwise convergent, as can be checked by estimating 
the ratio 



= ae r * a 



b n (b n +3H )t 



The left-hand side vanishes if 
3 



Then 



\br, 



-H 



1±\ 1 



2 m 
3#o 



(35) 



(36) 



Cf 



T{l + £-p) 



(£+l)T(£-p) 

- £, 1 - 2£ - 2 + p; -3H t) 
^(-£,l-2£ + p;-3H t) 
(£ - p)f(-l -£,l-2£-2+p; -3H t) 
(£ +!)*(-£, l-2£ + p;-3H t) 



n=l 



0. 



and one has local chaotic inflation (U(cj>) = const, the 
functions Oi vanish). The Friedmann equation deter- 
mines a set of algebraic relations between a n and b n co- 
efficients, leaving freedom to choose n* — 2 coefficients 



V. NONLOCAL SOLUTIONS: DIRECT 
CONSTRUCTION 



at large £, for finite Hq and t. The last step can be 
checked numerically or conveniently worked out by taking 
the dominant term n = £ in Eq. (|37p . For large £ and 
fixed p, the product can be approximated as 



2t-l 



U(p-k)~T(2£); 



k=0 



As said above, in general one has to specify one in- 
dependent Ansatz for each degree of freedom, (f>(t) and 
H(t). The procedure is inefficient and to find nontriv- 
ial solutions from this double prescription can be rather 
difficult. Here are a few examples. 



then T{2£ + 2)/[(l + £)T(2l)] ~ £. Hence j> is ill-defined 
and it is not a solution of the system. Note that, in 
order to discuss asymptotic solutions, the limit t — > oo 
must be taken after resumming the terms of the series, 
that is, once this is shown to converge. 
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B. Super-accelerating case (Ho > 1, q = 1), 4> — t p 

The l-th term in the series A reads 



□V = (-1)^6^- 



(38) 



h, n = (3HoY 



2n-l 



k=0 



-n in 



=0 d,,=0 \m=0 



xJJ(p-2i)^, 



(39) 



where the prefactor in front of t in the last line is a poly- 
nomial in (p — 2j), j = . . . n, of degree £ — n with all 
coefficients equal to f . In the hrst line, when n = the 
product is equal to 1. The coefficients bi tU increase with 
I and, again, is not well-defined for p ^ N. If p £ N, 
the series (|38[) terminates at n* = [p/2] (square brackets 
indicate the principal value of p: p = 2n* + I or p — 2n*) 
and <p is a hnite sum of powers of t, odd (even) if p is odd 
(even). In the case p = 1, for instance [4D solution for 
n = 2 in Eg. pgf]. 



t (X (j) . 



(40) 



C. Accelerating case (.Ho >1?9 — ~l))<^ = i p and 
<4 = lnt 



The first slow-roll parameter is constant, e = H , and 



□V = (-i) 



1^^-271)^-271-1 + 3^0) 



,n=0 



-21 



(41) 



If p G R \ N + , the coefficients of the series defining 
increase with l\ 



lim 



lim \d(p-2£) 



x(p-2£-l + 3H ) 



= +oo , 



for finite H and t. A possibility to avoid the convergence 
problem is that p is a positive even number, so that the 
series ends at £„ = p/2. Then 



p/2 
£=0 



21 



(42) 



is constant at early times, 6 ~ i p at late times, and the 
system is local (all the other functions Oi and U' are 



finite sums of powers of time). Another case where the 
series is finite is when 



P-1 + 3H = 2n , 



(43) 



for some natural n. 

The discussion easily extends to a more general Hubble 
parameter, q =/= ±1 in Eq. (|18p . We can immediately 
discard also the case </> = \nt since 6 would not converge: 



□V = (3H -1)(-1)' 



JJ 2(l + n)(2n + 3-3H ) 

."=0 

xt~ 2f , £>0. (44) 
The series is finite only when Hq = (2n + 3)/3, n E N. 

D. Preliminary assessment 

One may be tempted to conclude that either (a) the 
exact cosmological solutions described by Eq. (fP8|) can- 
not be used as a base for nonlocal functions, or (b) this 
Ansatz is viable only as long as the series defining <fi can 
be truncated at a finite order, in which case the system 
is trivially local at late times. 

As general statements, both (a) and (b) are proven to 
be wrong (also on Minkowski [5J). Although the defini- 
tion Eq. §3§ of the exponential kinetic operator is unique 
and unambiguous, what is not defined is the nonlocal 
solution expressed as an infinite series of powers of the 
d'Alembertian. We shall see in Sec. IVII Al that <j> defined 
in terms of the Green functions of the □ operator not 
only reproduces the special cases of finite series in t, but 
is well-defined even when the above approach fails. 

The reason is the following. We shall write an inte- 
gral representation for <fi which allows us to sum the ex- 
ponential series consistently. This solution satisfies by 
construction a diffusion equation which renders all equa- 
tions of motion local in time. Instead of having an infinite 
number of differential equations, one reduces the prob- 
lem to a set of algebraic equations in the parameters of 
the system. 



VI. NONLOCAL SOLUTIONS: INTEGRAL 
REPRESENTATION 

We propose a new recipe to construct nonlocal solu- 
tions. We interpret fixed value of an auxiliary 
evolution variable r, so that the scalar field <f>(r, t) is 
thought to live in 1 + 1 dimensions (in this discussion 
the 3 spatial dimensions play no role). The first step is 
to find a solution of the corresponding local system, that 
is with r* = 0. This can be done by solving Eq. © or 
(fTTI) with r* = 0, denoting the local homogeneous solu- 
tion 4>ioc(t) = 4>(r* = 0,t) and considering this as the 
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initial condition for a system that evolves in r through 
the diffusion equation 

ad r ^{r, t) = n^(r, t) , ^(0, t) = <f> loc {t) , (45) 

for some (unknown) parameter a. The solution of 
Eq. <g3J) is 

tf(r,f) = e^ a 4> loc (t) ■ (46) 

The action of nonlocal operators of the type e 9rD now 
becomes local for systems that solve Eq. (|45|) . In fact, 

e gra ip(r,t) =e aqrd ^(r,t) = V>((1 + aq)r, t) , (47) 

and all the effect of the operator is in the shift of the 
auxiliary variable r. The function tb(s,t), s ^ r, satisfies 
a new diffusion equation, given by Eq. (|4"5|) with r — > s 
everywhere (including the differential operator). 

Since r and d r do not commute, we need an ordering 
prescription for the exponential in Eq. (|47|) . We adopt 
the most natural one, setting all the derivatives d r to the 
right of the powers of r, in such a way that the translation 
property Eq. (|4"T|) holds also in this case. In addition, 
this is the only prescription compatible with the diffusion 
equation (j4"S"|) . In fact, 

e ra ^j(r,t) = ^□VM)=£^e? r V(r,i) 

= ^((l + o)r,i)- (48) 

In view of the above discussion, we can look for solu- 
tions of the nonlocal equations of the type </>(r, t) = 
C^{r)ij){r,i), where C,/,(r) is a constant in time. Since 
(j)\ oc is a solution of the local system (r* = 0), then 

CV(r = 0) = 1, (49) 

and we can choose C$ — e^ r , for some unknown (3. This 
is equivalent to allowing for a linear term in the diffusion 
equation for <p(r, t) : 

ad r </>(r,t) = a/3<f>(r,t)+n</>(r,t), (50) 

and the shift property (|4"T|) is unaltered. Equation (TTT1) 
becomes 

(□-™ 2 )Vj((7-a)r,t) =e~ 0r U'[e 0r ijj{{ 1 +a)r,t)], (51) 

where we started from a configuration (f>("fr, t), absorbed 
a factor 7 in /3, and used 

0( 7 r, t) = e /3r V((7 + oc)r, t) = e~ a/5 >(( 7 + a)r, t) . (52) 

The parameter 7 is not constrained dynamically like a 
and /3; it is introduced only to avoid any ambiguity in 
the present notation. One can fix it, for instance, so that 
7 = 1+a and cj>{r 7^ 0, t) on the left-hand side of Eq. (fijTj) . 

In general, the normalization constant Cjj in U will 
depend on r, so there is also the freedom to choose a 



function Cjj(r) such that 0/(0) matches with the nor- 
malization of the local potential. 

It should be stressed again that Eq. (f5Tj) is local in t, 
at variance with other methods developed in the liter- 
ature and based on a finite truncation of the exponen- 
tial defining the nonlocal operator [3^, HO, S2|. This 
property is very important, as it permits to study the 
nonlocal equations of motion in limited time intervals. 
Moreover, to understand the asymptotic (in time) dy- 
namical behaviour it is mandatory to look for solutions 
that, although approximate, sum the whole exponential 
series, rather than to find solutions at any fixed trunca- 
tion level D £ of the exponential operator. For instance, 
in Minkowski spacetime a fixed truncation of the spike 
solution Q gives rise to spurious wild oscillations near 
the origin, which are absent in the full representation. 

Finally, using Eqs. (|50|) and (|5Tj) , the functions 0\ and 
02 become (here 7 = 1 + a) 

O x = are 2rl3 f ds d Cr ^((r, t) {ad ir - m 2 )^(£r, t) , 

(53) 

° 2 = ,-2/n ^ i Ts I ^V>(Cr, t)(a%-m 2 ) 
V 1 ((1 + ot)r, t) Jo 

x far, t) , (54) 

where £ = 1 + a(2 — s) and £ = 1 + as. 

The calculation of (f>(r, t) proceeds as follows. 

(i) Find the eigenstates of the d'Alembertian operator. 
First, we solve the equation 

nG t ( P ,t)= ^G^t), (55) 

that is, 

Gi + 3HG % + fi 2 G 2 = , (56) 

where /i G C is either real or pure imaginary in or- 
der for □ to be a self- adjoint operator. The eigen- 
value equation is second-order and hence i = 1, 2. 

(ii) Write the nonlocal solution as an expansion in the 
basis of eigenstates ofd. The problem then is to 
find a function fi(/i) for each Green function Gi 
such that the local solution 4>(0, t) can be expressed 
as an integral transform, 

0(0, t) = J duldGxfat) h(p) 

+C 2 G 2 GM)/ 2 (M)], (57) 
while the localized field is a Gabor transform 5 

i>(r,t) = fdfie^/^dGxifi^Mfx) 

+C 2 G 2 (n,t)f 2 (»)}, (58) 



5 Strictly speaking, a Gabor transform is realized when the Gi 
are combinations of phases, which is the case of the Minkowski 
CSFT tachyon. 
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for some constants Cj, a and j3. These integrals 
become series if fi varies on a discrete set. The ad- 
vantage of the integral representation (|58|) is that 
the properties given by Eqs. (|45[) and (jT7|) are man- 
ifest. In other words, any function tp which can 
be represented as a convergent integral of the form 
obeys the diffusion equation. 



(hi) Verify that there exist some Ci, Cu, ct, and (3 such 
that Eq. (|58|) is a solution of the equations of mo- 
tion for r = rv This can be done in several ways, 
either by variational methods, or by expanding the 
equations in power series of r» or t and imposing 
the coefficients of the expansions to solve the equa- 
tions. 

The two approaches discussed in Sees. IIVI and IVl (covari- 
ant Ansatz and direct construction) are special subcases 
of this analysis. 

The heat kernel method of [B[ would seem more di- 
rect than the one based on the Green functions of 
the d'Alembertian operator. However, on curved back- 
grounds, the problem of finding a kernel satisfying a heat 
equation with covariant dAlembertian is much less triv- 
ial than in the Minkowski case. While on Minkowski 
spacetime such a kernel can be immediately guessed, on 
a FRW spacetime with a specific Hubble parameter one 
cannot proceed blindfolded and a systematic treatment is 
mandatory. An expansion of the tachyon field in terms of 
the □ eigenfunctions still yields the desired result with- 
out extra computational effort. 

We shall call local the system of equations solved by 
4>ioc(t) (r — r* = everywhere), while the (1 + l)-system 
solved by some 0(r, t) will be referred to as localized. In 
phenomenological models, there is freedom in the choice 
of the nonlocal potential in a way which does not affect 
the local equations of motion. Equation ([9]) can be writ- 
ten as 



U4> = h{4>)-h{4>) + U'{4>), 



(59) 



of motion depend on a finite number of slow-roll param- 
eters, and the slow-roll approximation is legitimate. The 
original nonlocal system contained infinitely many slow- 
roll parameters whose contribution could not be tested to 
be convergent; consequently, the slow- roll limit required 
much stronger assumptions out of analytic control [4l| . 

On a de Sitter background, the two independent eigen- 
functions are exponentials, Gj.,2 = e ±bnt , where b n is 
given by Eq. (|36[) with to replaced by /i. Since exact 
and nontrivial cosmological solutions require a varying 
Hubble parameter, we move directly to these scenarios. 



VII. NONLOCAL BRANEWORLD SYSTEM 
WITH H = Hot" 1 



In the power-law case, H = Hot 1 , the independent 
solutions of the eigenvalue equation (J5SJ) are ([Za], for- 
mula 8.491.6) 

G 1 (ii,t)=t"J v {fit), G2(»,t)=t"N„(tit), (61) 



(62) 



where J„ and N v are the Bessel functions of the first and 
second kind. 

Taking C'2 = 0, for /1 real and positive we can write 
the scalar field as 



r+00 

0(0, t)= / dnfi(n)Gi(ii,t)=1?F v (t) ; 
Jo 



(63) 



where 



Fu(t) 



(64) 



where h is an arbitrary function, while Eq. (HH) becomes is thc Hankel transform [77], whose inverse is 



- 2r ' ui = e- r * a [h(e- r * a j>) ~ h{4>)] + U'. (60) 



When h = and the nonlocal potential is simply the local 
one with the substitutions <f> — > <j), Uq — ► Cu(r), we will 
call the system described by Eq. (|59")l the minimal nonlo- 
cal system associated with the local one. The string case 
is nonminimal, h(x) = m 2 x. With this exception, there 
is no criterion for selecting h, and for the moment we 
will set it to vanish for simplicity (to 2 = 0). A nontrivial 
choice for h may become relevant when looking for exact 
solutions. Another kind of extra term in the equations 
of motion can be of the form r n h(4>, (j)), n > 0, which 
vanishes when r — > 0. 

The equations of motion of the localized system are lo- 
cal in time. This fact has two crucial advantages. First, 
there is only one momentum conjugate to the scalar field 
and quantization is well-defined. Second, the equations 



r+00 

flip) = M / dttF v {t)J v {,jLt) 
Jo 

r+00 

= fi dtt 1 ~ v cl)(0,t)J v (fj,t), 
Jo 



(65) 



for v > -1/2. 

In the braneworld (BW) case, using formula 6.561.14 
in M, 



/<Bw V)=/l(M)= ,„ni±^ r ^, ,„, 



r(i/-0/4) V2 



valid for -4 < 6 < —3H and v - 9/4 ^ —n, neff. 

The integral (|5"8)) can be done exactly when is 
Eq. (|66p. Using formula 6.631.1 of [7y|, the function 
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ipi(r,t) constructed from G\ is 

VMM) - * 2>-»"2I±^ 



r(i/-0/4) 



x / dfi^-^^e^ /a J v {nt) 



r(l + 0/4) / 4r 



9/4 / a* 2 ^ 



r(i/ + 1) V a 7 V 4r 
ai 2 



(67) 



where $ is Rummer's confluent hypergeometric function 
of the first kind, whose properties are summarized in the 
Appendix. The reader unfamiliar with it can just notice 
for the moment that it is defined as the series Eq. (|A.2|) . 
The integration domain is the region r/a < and 2v — 
6/2 = 1 — 3Hq — 9/2 > 0, but now one can analytically 
continue Eq. (|67|) for every Hq, r, and t (the typical values 
6 = ±1 do not constitute a problem). In order for the 
solution with r/a > to be real, one must compensate 
the phase e - m \ v -Q/v with the coefficient C\ in front of 
it. 

It is convenient to define 



so that 



ipi(r,z) = 4>o 



~Ar~ ' 



r(l + 0/4) ( Ar 
I>+1) 



(68) 



9/4 



— (-*)" 

a 



x9[v-~;v+l;z 



(69) 



Equation (f6"9")l satisfies, by construction, the diffusion 
equation (|45|) (see the Appendix). Notice that, for any 
and z, given < n, I S N: 

(1) v — 0/4 7^ — n and v + 1 ^ —I: tpi is a convergent 
series in z. If v + 1 = — i, the series Eq. (|A.2I) starts 
at fc = / + 1. 

(2) v — 0/4 = —n and i/ + 1 / — I: ipi is a polynomial 
of degree 0/4 in z. This condition is given precisely 
by Eq. (|4*5)) and it corresponds to one of the cases 
where the series representation of the e a operator 
is finite and well-defined, in the sense that the lo- 
cal solution lives in the operator's domain. This is 
an example where the analytic continuation to the 
region v — 0/4 < is explicitly required. 

(3) —0/4 = l — n+l, n > I, v noninteger: ip\ is an 
infinite hypergeometric series. 



(4) v - 0/4 = -n, v + l = -l, n> I: The series (|A~2j) 
has only the terms k = I + 1, . . . , n, with well- 
defined coefficients being T(k — n)/T(—n) finite. 
Again tpx is a polynomial of degree 6 / A = n — I — 1 
in z. This corresponds to p = 2,4,6,... in the 
notation of Sec. W\ 



(5) v - 0/4 = —n, v + 1 = -I, n < I: The Kummer 
function becomes an exponential times a polyno- 
mial of degree I — n = —1 — 0/4 in z, while the 
factor T(n — l)/T(—l) in front of it is finite. This 
case is unphysical (0 < —4) and corresponds to 
p = —2 (n = I) and p = —4, —6, —8, ... (n < I) in 
the notation of Sec. |Vl 

In the RS and GB braneworld cases, the last three in- 
stances are excluded. 

Some asymptotic behaviours of Eq. (|69|) : 

iz-i0(f->0orr-> oo): 

The function diverges at the origin of time if v < 
(H > 1/3). 

• z — ► oo (t — > oo or r — ► 0): Making use of the 
large z expression (|A.6[) . one has that, for z — > +oo 
(r/a > 0), 

_ +oc r(l + 0/4) / Ar\ 6/A 
Mr > z) ~ 0O I>-0/4) l"aj 

x(-I) ly e z z l/+1+e/4 , (71) 
while, for z —* —oo (r/a < 0), 



tpi(r, z) 



a 1 



•„ \ 1 



(72) 



A second, independent function which obeys the dif- 
fusion equation can be found directly from Eq. (j69]) . in- 
stead of defining and calculating an invertible transform 
with the eigenstate t v N„. 6 From what discussed in the 
Appendix, one can write rp2 as 



ip 2 (r,z) = 0o 



At 



6/4 



1 — V, Z 



(74) 



6 We have done it as a consistency check. Since G2(/J, t) = 
t v Nv(tit), one can write 0(0, t) as the TV-transform ( [78H . formula 
9.1) 

r+OC 

0(0, t) = 0o< fl/2 = / d^f 2 (ii)G 2 (iJ-,t), (73) 
Jo 

where ( 1791 . formula 11.1) 
/■+oo 

/ 2 ( M ) = M / dt 0(0, i)* 1 ""^ (/it) = /i(/x)tan(7T0/4) , 
Jo 

and H„ is the Struve function. The Gabor transform with the 
Neumann function N v turns out to be 

r + oo 2 

V>(M) = / dfie^ / Q / 2 (M)G 2 (M,t) 
Jo 

tan(7r0/4) , . , f tan(7r6»/4) 1 , , , 

tan(7r^) |^ tan[7r(^ — 0/4)J J 

where ip3 is defined below. 
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where $ is Rummer's confluent hypergeometric function 
of the second kind. An eventual phase is reabsorbed into 
C2. For any r and z, given < n, I G N: 

(1) v — 9/4 ^ —n and i/ + 1 ^ —I: ip 2 is a convergent 
series in z. 

(2) f — 9/4 = —n and i/ + 1 ^ — I: i/j 2 is a polynomial 
of degree 9/4 in z. 

(3) —0/4 = Z — n + 1, n > I, v noninteger: -0 2 is a 
polynomial of degree 9/4 in z. 

(4) -9/4 = l- n + l, v + 1 = —I, n> I: ip 2 collapses 
into ipi. 

(5) -9/4 = I - n + 1, v + 1 = -I, n < I; ip 2 is a 
polinomyal in z. 



Asymptotically: 



0/4 



I» 



r(-) 2 , 



r(i/-0/4) T(-9/4) 



i>2 4>ve™ e l 2 t e l 2 , 



02 



a e/2 . 



(75) 
(76) 
(77) 



According to the value of v, one should add next-to- 
leading order terms in Eq. (|75[) to get a consistent ex- 
pansion in z. At large z (small r), one recovers the local 
solution for a suitable coefficient C 2 depending on the 
sign of z (for general values of its first two arguments, 
^ has a branch point at r = 0). Modulo a phase, the 
local solution oc t 6 ! 2 is manifestly the initial condition 
at r = Q of tpi and ip 2 . 



Note that these properties hold for any r, including the 
analytic continuation into the r > region. Therefore, it 
is consistent to impose 



Ci = 0, 



(79) 



and take only one solution. 

In case 1 (—9/4 ^ l — n+1 and v+1 ^ —I), for generic 
values of 9 and v, the series built with terms as in Eq. (|4"Tj) 
is not convergent, while both 0>i and ip 2 are well-defined 
functions of time. When z < 0, ip 2 is complex- valued for 
general v. To get a real solution, there are two options: 
to set C 2 — or 



Ci _ ac i^-e/4) Sin(-7rfl/4) 



(80) 



It is worth recording this linear combination as some co- 
efficient C3 times 



ip3(r,z) = 0o 



T(l -I/ + 0/4) f 4r 

r(i-^ 



8/4 



with asymptotic behaviour 



2 _o , r(l-!/ + fl/4) / 4r\ e/4 

V3 ~ <po — — - — ( — r ) 



r(i-iz) 



Q 



(82) 



^3 
^3 



,-H-oo ^ r(l-!/ + fl/4) (_*r_\ 6,i e ^-t-8/4 



r(i-iz) 



(83) 
(84) 



A. Comparison with the series calculation 

Before looking for a solution of the nonlocal equations 
of motion, we compare the function 



ip = Clip! + C 2 l\) 2 , 



(78) 



found via the eigenfunction expansion method, with the 
brute- force series computed in Sec. IV Cl First, we check 
that it is always possible to choose a suitable linear com- 
bination of Eqs. dni]) and ([73]) reproducing the series 
when this is finite (that is, when ijj is a polynomial). 
Then, for cases 2 to 5 (a — 1 in what follows): 

(2) Both ipi and ijj 2 are proportional to the polynomial 
found via the series computation. 

(3) Fixing C\ = 0, ij) is proportional to the polynomial 
found via the series computation. 

(4) Same as case 2. 

(5) Same as case 3. 



Figs. [T] and O show the function i/j(r, t) as superposition of 
eigenstates of the d'Alembertian and as truncated series. 
The dashed line is the solution of the local system (r = 
0). The dotted curves are the series functions given by 
Eqs. © and ©, truncated at I = 2,5,8,11 (from left 
to right). The solid curves are linear combinations of ipi 
and tp2 as explained in the caption; these are the main 
result of this part, to be compared with the truncated 
series and the local case. To summarize: 

(i) Case 1 (v - 0/4,1/ £ Ir): For r/a > 0, only 2 
reproduces the series computation at large t, while 
for r/a < only a real linear combination of ipi and 
ip3 is allowed. Note that the boundary condition 
given by the local solution of the local system is 
automatically satisfied. 

(ii) Cases 2 to 5 (finite series): Only ip 2 , for any r/a, 
is required. In case 2, being the physical one, ipi oc 

02- 

(iii) In all cases, the heat equation is fulfilled. 



11 













/ 








1 


2 


3 4 5 6 
t 




-1 



\ \ 

\ \ 
- \ \ 
\ \ 

s.' 

















t 



FIG. 1: The function tf)(r,t) calculated via the series ex- 
pansion and the eigenstates method. Here, 9 = 1 (RS 
braneworld) and v = —3/2. The meaning of each dashed 
and dotted curve is explained in the text. In the upper panel 
(r/a = 1), the solid thin curve is if>i with Ci = e ^-e/4) ^ 
while the solid thick curve is ip2 with G% = 1. In the lower 
panel (r/a = —1), the solid thin curve is f/>i with Ci = 1, 
while the solid thick curve is ips with C3 = 1. 



FIG. 2: The function tf)(r,t) calculated via the series ex- 
pansion and the eigenstates method. Here, 9 = — 1 (GB 
braneworld) and v = —3/2. The meaning of each dashed 
and dotted curve is explained in the text. In the upper panel 
(r/a = 1), the solid thin curve is if>i with Ci = e i7r( - v - 0/i) , 
while the solid thick curve is -i/>2 with G% = 1. In the lower 
panel (r/a = —1), the solid thin curve is f/>i with Ci = 1, 
while the solid thick curve is 1^3 with C3 = 1. 



B. Solutions of the equations of motion 

The problem of finding solutions of the scalar nonlocal 
equation with the local solution as boundary condition 
gets simplified if one takes into account only one of the 
three functions ipi discussed above. Considering z and 
r as independent variables and using the diffusion equa- 
tion, the equation of motion for the scalar field is now 
first-order for solutions of the form tp(r, z) = r~ b u(z) for 
some b. Actually, 

Dip = -(rd r - zdM = --(b + zd z U . (85) 
r r 

Eqs. (TA~1~2|) , |A~13)) , an d (gU) yield 

^+(r, t) = -- Cu(r) e-^l e ((1 + 2a)r, t)} 1 ^ 9 , 
a 

(86) 

and 

1#(r, t) = ^ Cu(r) e~^/e ^ ((1 + 2a)r , t)]l -4/fl ; 

(87) 



while for -03 

^(r,t) = -(4i/-0)C c/ (r)e- 4/3 ''/ e 
a 

x& 3 ((l + 2a)r,t)] 1 - 4/e , (88) 

where the function Cu (r) is arbitrary except for the prop- 
erty Cu — > 1 as r — > 0, and the superscript + indicates 
that the first argument of the Kummer function in ipi is 
translated by +1. Case 1 with r/a > and case 2 for 
any r/a are described by Eq. (|87|) , while for case 1 with 
r/a < one can choose between Eqs. (|8"6"T) and ([55)) . Case 
1 with r/a < also admits a linear combination of ipi 
and ^2- 

Here one can clearly see the nature of the Cauchy prob- 
lem in nonlocal systems. These equations are purely al- 
gebraic in the scalar field. To solve the original equations 
of motion, one should find an infinite set of initial condi- 
tions 4>q, <po, (£>o, . . . . If 4> is analytic, this is equivalent 
to find the Taylor expansion of 4> around a point, and to 
know the initial conditions means to already know the 
solution [57J . In the above expressions, the set of initial 
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conditions has been reduced to a finite one given by the 
parameters fa, Hq, ... of the solution cj>, which is known 
in its functional form thanks to the heat equation. 

In the original equation of motion for <j) there appeared 
all the derivatives of the field, that is, the full tower of 
slow-roll parameter. While in local systems there are only 
a finite number of slow-roll parameters and all these can 
be set to be small, in the nonlocal case it is not guaran- 
teed a priori that the slow-roll tower as a whole is negligi- 
ble with respect to 0(1) terms; the results of Sec. IVII Al 
are a demonstration of this fact. In [4l[ an extra con- 
dition for this approximation to be true was given, but 
here we do not even need it. In fact, <j)(r, t) depends only 
on the first slow-roll parameter e = 3/(1 — 2v), and the 
standard slow-roll approximation can be safely applied. 
This is a nontrivial property which may justify a posteri- 
ori the truncation scheme sometimes adopted in the past 
[39T |40| . |42| in the case when the full nonlocal solution is 
known. 

Let us now find solutions of the dynamical equations. 
First, we consider Eq. (|57|) . An exact solution can be 
inferred by looking at the asymptotic behaviour of fa; 
expansions in z (that is, t or r) are always well-defined 
and there is no issue of convergence (although there may 
not be a global solution of the form <j>(r,t)). At large 
z, matching powers of z cancel out and at lowest order 
[M = in Eq. ([ATTT]) ] one gets 



Cu(r, 



_ e H3r/9 



(89) 



The scale is always fixed if Cu(r*) is given from the 
beginning or, as alternative, C\j defines the shape of the 
renormalization factor in front of ip(r,t). At next-to- 
lcading order (M = 1), the equation of motion is safistied 
for a = 2/(4v—9), while at third order (M = 2) one finds 
the conditions v = 1 (describing a contracting universe, 
which we will ignore) or 



(90) 



giving a case 2 exact solution of Eq. ([57)) in high-energy 
braneworld power-law cosmologies: 

fa{r 1 t) = U e/2 {l+^) ■ (91) 

In the Randall-Sundrum braneworld (6 — 1), v — 
-3/4 =*> H a = 1 - 6/6 = 5/6 < 1, which is out- 
side the accelerating regime. In the Gauss-Bonnet case, 
v = —5/4 Ho = 1 — 6/6 = 7/6 > 1 which gives very 
mild acceleration. This solution is steeper than the local 
one near the origin, as one can also see by considering 
the expansion around z = of fa, taking — v 3> 1 and 
noticing that powers of z do not match in the left- and 
right-hand side of Eq. (|87|) . The sign of r6 determines 
whether <f> — > ±00 as t — > 0; if one wants to stick to the 
r > case, the GB <j) is an increasing function which 
vanishes at t* = \/2r. 



The value a = —1/2 may seem problematic because 
the local function with r = in the right-hand side of 
Eq. (|87| would be equated to a nonlocal part with r/0. 
This is only a matter of notation. In our formalism, r can 
be translated arbitrarily via e a operators, and the equa- 
tion of motion remains nonlocal regardless what transla- 
tion is chosen to represent the r dependence. The special 
value a = —1/2 defines a "coordinate system" where 

4> cx 4>\ oc . 

The field <j> in the Friedmann equation is e@ r fa{r /2, t). 



When a = -1/2, Eqs. (53]) and (H]) become (7 



0) 



Ox = ~re 2 ^ r dsd sr fasr/2,t)d sr fa(l-s/2)r,t), 
Jo 

(92) 



On 



2r 



fa(r/2, t) Jo 



ds fasr/2, t) d sr ip((l - s/2)r, t) . 



(93) 



Equation (|91|) is linear in r and the computation of these 
quantities is straightforward: 



0i oct* 



t 4 



[t 2 + (6- 4)r]' 



1 



(94) 



The Friedmann equation is / F Rw(i) 



t v 



[O(l) 



0(t~ 2 ) + 0(t~ A )] and is not solved exactly, since the Oi 
do not match with the the rest of the equation, ~ t e ~ 2 . 
It is possible to choose such that the O(l) term in the 
Friedmann equation vanishes, but the other terms cannot 
be cancelled. At large t the Oi are subdominant and the 
usual local cosmology is consistently recovered. It is con- 
firmed [1^, |4l| that, in principle, the sign of the kinetic 
term of <f> can change during the evolution of the universe 
(at early times, O2 ~ — 1, while at late times O2 ~ 0). 
For instance, if r < 0, sgn(l — O2) > always; if r > 0, 
sgn(l— O2) < (and actually blows up) for a finite period 
during the evolution. However, this qualitative result is 
not robust because, at times when sgn(l — O2) < 0, the 
Friedmann equation is not solved with a satisfactory level 
of accuracy, \f FRW (t)/(H 2 - e + 2 p)\ > 10%. 

Equation (l9~lj) is also the only asymptotic solution of 
the scalar equation at early times. 7 Since it is not an 
early-time solution of the Friedmann equation, there are 
no nontrivial solutions of the nonlocal system at early 
times. 

To get solutions falling in case 1, one has to con- 
sider more general values of v. For instance, — v 3> 1 
corresponds to slowly rolling fields in an expanding uni- 
verse. We have checked for all fa and several values of the 



7 The reader can check this by expanding the equation of motion 
with a = —1/2 for t <§C 1. ipi collapses into the same function, 
while tp3 ~ const + 0(t 2 ). One can absorb the constant in the 
definition of ip3 but it is not possible to match the next-to-leading 
term in the right-hand side of the equation of motion. 
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parameters that none of equations of motion is exactly 
solved, and the only regime where they are solved approx- 
imately is at late times, that is, when <j){r, i) tends to the 
local solution t e l 2 . To be more precise, the equation of 
motion for the scalar field is of the form A = B, where 
A and B are noncomposite functions of time. One can 
choose a time f* above which \(A — B) / (A + B)\ < 10~ n 
for some given n. A way to get approximate (and still 
almost local) expanding solutions at early times is to con- 
sider v > and then perform a time reversal of the equa- 
tions. If v — 9/ A is chosen to be a positive integer, <f>(r, t) 
can be written in terms of exponentials, powers and in- 
complete gamma functions. 

We note that the higher the value of |f|, the higher the 
characteristic time t* . This important feature is compat- 
ible with the notion that it is difficult to realize inflation 
in nonlocal theories, with the possible exception of the 
special case, described in Sec. lIVl where D<f> oc <f> |44ll4Sj]. 
If analytic, the solution of the equations of motion is 
uniquely determined by the infinite set of initial condi- 
tions one should specify in the nonlocal Cauchy problem. 
As we have seen, the latter can be conveniently translated 
into a finite set of algebraic constraints, which fixes most 
or all of the parameters of the theory. The solution is 
very rigid and not all values of the parameters are al- 
lowed; it firmly holds in memory the initial conditions, 
and any mechanism (such as inflation) which tends to 
wipe them out is unlikely to be realized. 

The condition a = —1/2 simplifies the calculations be- 
cause in that case the potential term is a monomial, but 
we have verified that all these results can be generalized, 
with no appreciable changes, to other values of a. 



rectly and reads 
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(96) 



2I» ^ k\ 



(97) 



which is convergent with infinite radius of convergence. 
When v S Z _ , the series terminates at k = fc max = —v. 

Here we point out a trick which allows to get the same 
result with slightly less computational effort (in fact, 
the direct calculation relies on a similar shortcut). One 
should treat the 4D and braneworld cases separately be- 
cause, setting 9 = in Eq. (JMU, / (BW) | e=0 + / (4D) - This 
is expected, since the local braneworld solution <fi oc t e / 2 
is trivial for 9 — 0. However, the following relation be- 
tween the local solutions is valid: 



Au( Bw y/2] 



0=0 
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2^ 4D) 
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(98) 



where we have treated 
can verify that 



,(BW) 



as a fixed constant, and one 



VIII. NONLOCAL 4D SYSTEM WITH H = Hot' 1 

In four dimensions (</> = </>o hit), using formula 6.771 in 
[76|, the kernel realizing the integral transform Eq. ([63]) 



/ (4D) M 



h(p) 

4, (*)"[«,, 



(95) 



where ip = c^lnl^a;) is the digamma function and 
7 = —ip(l) ~ —0.577 is the Euler-Mascheroni constant. 
The integral (|65p is done in the region of the parameter 
space n > 0, v > 1/2, implying H < 0. One can ei- 
ther analytically continue Eq. ([93)1 to the region v < 1/2 
or, if interested in post-big bang cosmology, make a time 
reversal to get a superaccelerating solution. 

The resulting nonlocal function can be calculated di- 



/ 



,» (/i)= 2^ D; 5/( BW )( M ) 



,(BW) 
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09 



(99) 



)=0 



Then one can compute <fi(r, t) in the simpler case of the 
braneworld and then get the nonlocal field in four di- 
mensions by just a derivation. We need formula (|A.14|) . 
yielding 



d e <$>(v 



9/A-\ + V - z)\ g=0 

k) — ip(y) z k 



+oo . , 



E 

k=l 



k 



k\ 



after 



Taking the derivative with respect to 9 of Eq. (|69|) 
the rescaling of Eq. (|99|) one gets exactly Eq. (|96l) . 

The asymptotic behaviours of tpi can be calculated di- 
rectly or by taking the ^-derivative of Eqs. ([70)1 - (|T2"|) . In 
particular, 



Vi 



V(i) 



1 



In 



2T(v + l) 
while in the limit z — > +oo, ip\ cx e z z 



1 (const 



(100) 
■In z). 
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Although one can define an antitransform via a func- 
tion fi(ij) oc fi"^ 1 , the integral transform Eq. (j?3"j) for 
G2 is singular, as the parameter range lies outside the 
definition domain of the transform. 

Nevertheless, in the braneworld case this transform ex- 
ists, and the resulting functions ip2,3 are well-behaved. 
Another 4D nonlocal function can be found by deriving 
one of the other independent braneworld functions with 
respect to 8. For instance 

d e $(-e/4; 1-u; z)\ g=Q 

r(i - v) r(fc) z k 

~ 4 T(l - i/ + A;) IT 

2 F 2 (1, 1; 2, 2 - u; z) , 



4(1,-1) 



and one gets 



1 /4r 
l \ a 



2 F 2 (1, 1; 2,2-iz; z) 



(101) 



where 2-F2 is a generalized hypergeometric function de- 
fined in Eq. (|A.16p . V>3 is well-defined only when v < 1. 
Asymptotically, 



1p3 



z->0 90 



1 /4r 
^ (1-1/ ) + -In - 

2 \ a 



0o In t , 



(102) 



(103) 



while in the limit z — * +00, 1P3 oc e z z I/_1 (const — lnz). 
It is checked that Eqs. (|96|) and (|101[) obey the diffusion 
equation (1431) . 

A particular linear combination which can be derived 
from the braneworld function ip2 is 



1p2 = 1p3 ~ y 



+ r(-Z/)(-z) I/ $(z/; l + i/; z) 



(104) 

This function differs from tp 3 by a term which vanishes 
at late times, and generalizes the constant field shift 
parametrizing the class of exact local solutions. 



A. Comparison with the series calculation 
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FIG. 3: The 4D function ij}(r,t) calculated via the series ex- 
pansion and the eigenstates method. The dashed curve is the 
local solution cj> = lnt, while the dotted curves are the series- 
type functions at different level truncations. The solid curves 
are (increasing thickness) ipi, ip2, and i/>3, respectively. Here, 
v = -3.9, r = 0.26, a = -1. 



B. Solutions of the equations of motion 



The local potential is U{4>) = L/ e~ 2, ^ Q , where U = 
—v<j)Q and cf>o — ^/2Hq/ n 2 . Under the substitutions <ft — > 
(j) and Uq — > UqCu{t), the nonlocal equation of motion 
for the scalar is 



□V>i(r, t) = 2v e~ tir Cu(r) exp[-2e /3 > l ((l + 2a)r, t)} , 

(105) 

where we have cancelled the <fio factors in the ipi's. The 
function Cu(r) is arbitrary except for the property Cjj — > 
1 as r — > 0. To get the local solution as an asymptotic 
solution at late times, it must be Cu{r) — e^ r ■ 

The left-hand side of the equation of motion can be 
worked out explicitly. It is uninstructively complicated 
for ipi, while for ip2 and ip3 



(106) 

Oh = f {1 + eVT(l -v)- r(l - v, z)}} , (107) 
2r 



Figure [3] shows the functions if>i(r,t) together with the 
divergent truncated series calculated before. The dashed 
line is the solution of the local system (r = 0). The 
dotted curves are the series function given by Eqs. ([3]), 
(US!), and 05), truncated at £ = 16,18,20,22 (from left 
to right). The solid curves are ^1,2,3 as explained in the 
caption. 

Again, it is possible to choose the parameter v (H ) 
in a region where the series function diverges while the 
functions ipi (r, t) calculated via the diffusion equation are 
well-defined. 



where E v {z) = ds s u e sz is the exponential integral 
function and T(l — v, z) = ds s~ v e~ s is the incom- 
plete gamma function. For r/a < 0, the first expression 
is real, while for r/a > the definition (|104|) needs an 
extra (— l) v factor in the last term, so that the term in 
square brackets in Eq. (|106[) vanishes. 

In order to find asymptotic solutions for general z/, one 
can choose a = —1/2 so that the potential becomes a 
power oft (t~ 2e ) which can be matched with the left- 
hand side of the equation. Assuming v < 0, Eq. (|105|) 
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for ipi is satisfied at lowest order in t <C 1 when 



e? r = -v. 



(108) 

° U ~ (8r)^r(l + ,) ■ (1 ° 9) 

The sign of Cjj depends on the choice of v. For 
Eq. (fTU5|) still holds together with 



Ci 



r(i-^) 
( 8r )i+^ 



(110) 



In an accelerating regime with r > 0, Cr/ < 0. As re- 
gards ip3 at early times, at next-to-lowest order -03 ~ 
— l/(4r) + 0(t 2 ). One can absorb the constant in the 
definition of ip^, and eventually match the next-to-leading 
term. The latter step can be done after a generalization 
of the renormalization factor e@ r to a function (r) such 
that C^(r) — > 1 as r — * and C^,(r*) < 0. 

The Friedmann equation is not solved at early times by 
these choices of the parameters, and as in the braneworld 
case the only acceptable asymptotic solutions are those 
at late times. 



IX. DISCUSSION 

A. Summary of the main results 

We have presented a method which allows us to explic- 
itly work out nonlocal fields of nonlocal systems which 
do not necessarily admit a representation as a series of 
□ operators. This is the main difference relative to other 
proposals: it is not mandatory to truncate a series-type 
solution to calculate its coefficients, because the solution 
itself is already resummed. The simple case of power- 
law cosmology has been studied both in a pure four- 
dimensional scenario and on a brane with modified high- 
energy Friedmann equation. In all cases the found solu- 
tions are indeed approximate and close to the local be- 
haviour. 



B. Comparison with other studies 

We want to compare our results with those recently 
appearing in literature. We already went through the 
subtleties of a truncation procedure [H, 52, so we 
comment on the remaining points of view proposed on 
this stimulating subject. 

In Ref. [43|, [4f|, the operator e D is represented as a 
Weierstrass product which allows one to recast the non- 
local tachyon Lagrangian as that of an infinite number 
of decoupled scalars ipn which are eigenfunctions (with 
eigenvalue /i 2 ) of the d'Alembertian. The latter is not a 
self-adjoint operator acting on these scalars, and the re- 
sulting total Hamiltonian has a continuous real spectrum 



which cannot describe particle modes. 8 

However, the □ operator is assumed not to be self- 
adjoint because the authors look for kink-like solutions 
mediating from the unstable open-string vacuum to the 
stable closed-string vacuum at t — +00. Then, the char- 
acteristic equation for the differential operator cannot be 
solved for finite (real) eigenvalues /it 2 in Eq. (|55|) . As a 
consequence, solutions of this kind necessarily oscillate. 
Note that if one extends the approximate solutions of 
[42l [45l | also to negative times, these may have wild os- 
cillations. 

On phenomenological grounds, the author of [Ifl con- 
siders a class of Lagrangians of the form C = 4>F(E\)(f> — 
a<f> p+1 , where F is an unspecified complex operator. A 
notorious t oy model falling into this class is the p-adic 
string [79|, Ho, [H|. The analysis assumes that F(z) is 
analytic on a certain domain including the point z = 0, 
so that F admits Taylor expansion around the origin. 
These requirements translate into imposing that the op- 
erator F(\3) can be written as a power series in □ when 
acting on the physical states <j). The field <f> is written 
as a superposition of the scalars ip n , which have conse- 
quently the same role of the Green functions G(/i, t) of 
Eq. ([55]) . Here we have considered only the case /i 2 6 E, 
but it would be interesting to apply our formalism also 
to cases with Im(/i 2 ) 7^ and confront the results with 
the perturbative solutions of [45j. 

In [46| , an action with a potential of the form (|T2]) with 
a quadratic nonlocal term is solved exactly on Minkowski 
space (via a Weierstrass product formulation) and its so- 
lutions of type J2 n C n e^ nt are classified. On an FRW 
background, solutions are found under three different ap- 
proximations: a truncated expansion (up to some n max ) 
in terms of the local scalars ip n which are also equipped 
with ad hoc potentials; a truncation of the Weierstrass 
product; and a truncation of the Taylor series defining 
the kinetic nonlocal operator. Either method is effective 
only for computing finite (eventually large-order) series- 
type solutions; the order of the truncation determines 
the number of spurious solutions of the equations of mo- 
tion (see also [4l| for a discussion on this point). Also, 
it is not clear whether the same treatment can lead to 
exact or well-defined approximated solutions in the case 
of more realistic potentials like those of effective SFT. 

The p-adic string in a cosmological context is studied 
in (44l [49l | , where a relaxation of the ?y problem is pro- 
posed (see also |48j|). Again, the considered solutions are 
finite series and it is conjectured that resummation of 
the kinetic operator is possible at arbitrary order. Here 
we gave negative examples with a different Ansatz for </>, 
so one cannot disprove this statement. Even if resum- 



8 See [47l for a first step in the treatment of the quantization prob- 
lem. Note that the 1 + 1 Hamiltonian formalism avoids all diffi- 
culties associated with an infinite number of degrees of freedom; 
namely, there is only one local, physical field (j>(r, t) to be quan- 
tized via a finite set of canonical commutation relations. 



16 



mation were not possible, one could not conclude that 
this type of solutions is not viable; according to our per- 
spective, it would not belong to the spectrum of the ki- 
netic operator written as a series, and another handling 
of the latter would be required. In [44|, suppression of 
the second-derivative term in the d'Alembertian relative 
to the Hubble friction is assumed, and e D w e -3H d t _ 
This approximation is valid only by virtue of the special 
kind of solutions considered, while it is not in general. 



C. Problems and open issues 

The main virtue of the integral representation of non- 
local fields is to show that nonlocal theories are tractable 
in a quantitative way. However, the analytic and numeri- 
cal inspections of the dynamical equations for the chosen 
example (a = t H °) are still incomplete and the fact that 
one or both equations of motion are not exactly solved 
may be subject to different interpretations. 

(i) Extra (decoupled) particle fields might be in- 
voked and their evolution equations integrated 
into the system, as for level corrections to the 
CSFT tachyon potential. However, such eventual- 
ity would not escape the problem of finding single- 
field nonlocal solutions given a local system as ini- 
tial condition. 

(ii) The failure of the equations of motion might in- 
dicate that the local solution cannot be an initial 
condition for the nonlocal system. Other initial 
conditions might be concocted but this case would 
undermine one of the basis of the proposed method; 
that is, the possibility to attack nonlocality thanks 
to the systematic knowledge of the physics at r = 0. 

(iii) Third, we have not considered modifications of the 
minimal system with "evanescent" terms of the 
form h(4>) — h((f>) or r n h((j>, 0), n > 0. These can be 
tuned so that to adjust the asymptotic behaviour at 
early times and may provide a necessary ingredient 
for the solution of the dynamical equations. 

(iv) A fourth possibility is that the nonlocal cosmolog- 
ical system defined under the substitution U(cj>) — > 
U(cj)) in the local equations has no (analytic) solu- 
tion except the local one. We do not know whether 
this "no-go" conjecture is a general characteristic of 
nonlocal FRW systems or only of those with power- 
law scale factor. Counter-examples on a noncosmo- 
logical (in particular, Minkowski) background are 
already known 0, || . 

The space iSdc of functions obeying the diffusion equa- 
tion includes states found via the direct computation 
method in the cases where the series is finite (i.e., for 
a particular choice of the parameters); we denote this 
subspace as <S*. Therefore: (A) either the space iS p h ys 
of physical states (that is, those satisfying the equations 



of motion) is a subset of the direct approach for finite 
series (S P hys C «S»), or (B) the solutions obey the diffu- 
sion equation but are infinite series (<S P hys C Sdc \ <S*) , or 
(C) they lay outside Sdc- While (C) is beyond the scope 
of this paper, we did not yet found any conclusive result 
about (A) and (B). In fact, the particular examples we 
considered (power-law universes) seem not to have any 
nonlocal nontrivial solution, but there is no evidence that 
such a negative situation is a general feature of nonlocal- 
ity on FRW. On the other hand, solutions of the string 
equations of motion do belong to case (B) [J, Q . 

There are several open problems which should be part 
of the future research agenda. Among them: 

(i) Models with other choices for the Hubble param- 
eter and the local potential should be considered. 
Superaccelerating examples [Eq. {23])] and the case 
of the CSFT tachyon are of particular relevance 

(ii) The dynamical system can be studied also within 
the Hamilton-Jacobi formalism, where the coordi- 
nate t is traded for the scalar field as a measure of 
time. In particular, the Friedmann equation fT^ 
can be simplified into [4l[ 



e(2-0) 
6 



H 2 ~ e 

— — + Oi-V = 0, 



4 



(in) 



where the dependence on Oi and cj> has disappeared 
and <j) is regarded as the natural clock of the sys- 
tem. A preliminary check shows that our numerical 
results are not modified in this framework. 

(iii) The classical (i.e., phase-space) stability of the non- 
local solutions is to be assessed. 

(iv) Cosmological linear perturbations would be impor- 
tant to describe inflationary observables as well as 
semiclassical stability of nonlocal models. Even if 
the Minkowskian theory is ghost-free, this is not 
guaranteed on a curved background. However, 
we expect that small-scale perturbations are un- 
changed in nonlocal scenarios. Consider the pertur- 
bation equation Eq. (4.41) of 4l| for a test field (no 
metric backreaction) . When the comoving spatial 
momentum of the perturbation is large, |k| — > oo, 
the effective mass term is subdominat and the har- 
monic oscillator equation of standard inflation is 
recovered. 

(v) We have not said much about the singularity 
problem. Are there nonlocal dynamics which are 
smooth at the big bang? 

(vi) Since the scalar field behaves as a perfect fluid, 
Eq. (fTB)) relies upon the same assumptions advo- 
cated in the local braneworld. The first is that 
there is a confinement mechanism such that mat- 
ter lives on the brane only, while gravitons are free 
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to propagate in the empty bulk. This is guaranteed 
as long as the brane energy density p is small com- 
pared to the effective gravitational 5D mass cou- 
pling. Then the energy-momentum tensor is co- 
variantly conserved and the continuity equation, 
i.e., the scalar equation of motion, has no source 
term (correspondin g to the 05-component of the 
bulk tensor [H, H lH US HE M, 113 )• This 
does not imply that the contribution of the Weyl 
tensor is neglected (the second assumption): given 
a standard continuity equation on the brane (empty 
bulk), the Friedmann equation still can get an ex- 
tra dark-radiation term. In the local braneworld, 
bulk physics mainly affects the small-scale or late- 
time (i.e. post-inflationary) cosmological struc- 
ture and can be consistently neglected during the 
early history of the universe. For local inflation, 
this is a highly nontrivial result which has been 
confirmed with several methods both analytically 
and numerically ( [5lT | is among the most recent in- 
stances). Intuitively, the Randall-Sundrum dark- 
radiation term, which is the simplest contribution 
of the Weyl tensor, scales as a~ 4 and is damped 
during the expansion of the universe. 

If the gravitational sector is local, these results do 
not depend on the matter content as long as in- 
flation is enforced. If it is not (a possibility in- 
dicated by our findings), only an investigation of 
a concrete model would set the issue. However, 
a dark-radiation term in the Friedmann equation 
would not affect the shape and parameters of the 
nonlocal solutions, because these are determined by 
the scalar equation. In principle it can modify the 
early-time behaviour, but in that region the func- 
tions we found cease to be solutions of the Fried- 
mann equation. We checked this even in the pres- 
ence of a simple Weyl contribution oc a~ 4 = t~ AH ° . 
In this case, moreover, the accuracy of the solution 
is considerably improved at intermediate times and 
the confidence value i* (Sec. I VII B[) is decreased. 

The validity of the simple configuration of a brane 
with nonlocal 4D matter in an empty 5D bulk de- 
pends on the details of the model and cannot be de- 
cided by heuristic considerations. One of the prob- 
lems is the implementation of nonlocality also in 
the gravitational action (see below). 

(vii) It would be interesting to quantify the relation be- 
tween the heat equation approach and the 1 + 1 
Hamiltonian formalism [12] . The main difference is 
that the diffusion equation is second order in time 
derivatives and first order in r-derivatives, while 
the auxiliary variable in the 1 + 1 formalism acts 
as a time translation (first-order equation). So it 
should be understood how to incorporate the diffu- 
sion equation as a covariant constraint in the La- 
grangian. 

(viii) We already mentioned the necessity to extend our 



construction to instantonic solutions. This may be 
related to another aspect constituting new terri- 
tory of research in cosmology: boundary string field 
theory (BSFT), for which solutions now have been 
constructed which are in deep relation with those 
of CSFT [5[. One can extend the calculation of 
the BSFT energy-momentum tensor of [901 ] to an 
FRW background and study the resulting Fried- 
mann equation. 



(ix) It should be underlined that to put a curved met- 
ric into the original effective action of CSFT is 
not as straightforward as it may seem. A grav- 
itational background is mediated by the 2-tensor 
modes of the closed-string perturbative spectrum 
(defined on the stable vacuum), while in the ini- 
tial configuration only open-string modes (includ- 
ing the tachyon) live on the unstable brane. The 
direct substitution t]^ — > g^ v into the open-string 
effective action is certainly an intuitive step help- 
ful in the description of cosmological scenarios, but 
further meditation on this issue according to the 
rules of string theory would be of benefit. 



(x) For consistency, one should consider the effect of 
nonlocal interactions also in the gravitational sec- 
tor and check how the dynamics discussed in this 
paper is modified. The nonperturbative string for- 
mulation of gravity is still unknown, and only phe- 
nomenological Ansdtze have been proposed so far 
[75l li)H . One possibility would be to adopt an ed- 
ucated guess about the form of the effective grav- 
itational action and extend the heat equation for- 
malism to its degrees of freedom. The task seems 
lengthy but at hand; of course, its results would 
be sensitive of the choice of the action. Another 
avenue to explore might be to truncate the non- 
local scalar part to finite order in a' and cou- 
ple it to gravity at the same order. However, it 
seems unlikely that a perturbative approximation 
would shed new insights in the problem. The rea- 
son is that, due to higher-order differential oper- 
ators, ghosts and other instabilities would appear 
in all sectors. These instabilities have been stud- 
ied at second order in a' (Gauss-Bonnet modulus 
gravity) for an FRW background in [92j (see also 
[H, 0, [sl, [96j])- One way to escape them might 
be to consider very realistic models with particular 
potentials and matter content. All of these ingre- 
dients should be dictated by string theory and un- 
der careful control throughout the evolution of the 
universe. The price to pay for abandoning the non- 
perturbative regime would be an extremely com- 
plicated model whose features would be accessible 
only numerically. 
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APPENDIX: HYPERGEOMETRIC FUNCTIONS 

We collect some properties of hypergeometric func- 
tions, widely used in Sees. IVlil and I Villi In the following 
and except in Eq. (|A.5[) , a and j3 are arbitrary constants 
which have nothing to do with the coefficients in the non- 
local equations of motion. 

The equation 



''dz 2 



fill 

(p-z) — -au = 0, (A.l) 



admits two solutions ([Za|, formulas 9.216.1-3), which are 
Rummer's confluent hypergeometric functions of the first 
and second kind: u\ — <J>(a; (3; z), u-2 = \t(a; /3; z). u% 
(also called iF\ or M) is defined as 



$(a; p; z) 



+oo , , u 

\ - (a)k z 



^ r(a + k) 



m 

r(a) ^ o r(/? + fc) k\ 



where (a)k = T(a + k)/T(a) is the Pochhammer symbol. 
Rummer's function of the second kind \& (also called U) 
is: 



sin 7r/3 



3(a; p\ 



r(l + a-/3)r(/3) 
$(l + a-/3; 2-/3; z) 



r(a)r(2-/3) 



(A.3) 



which satisfies the identity 
#(a; /3; z) = z 1-/J #(a! 



/3+1; 2-/3; z) . (A.4) 



To link this discussion with the physical system under 
study, we notice that the diffusion equation can be writ- 
ten as an operator acting on the two independent vari- 
ables z [Eq. flBSJ)] and r: 



-{ad r - n)=zdl + {l-v- z)d z + rd r 



(A.5) 



which reproduces Eq. (|A.1[) for u(r, z) = r a u{z). The 
solutions jM]) {a = v - 9/4, = 1 + u) and |j72J| (a = 
—0/4, /3 = 1 — v) are of this form. 

The properties of Rummer's functions when a and/or 
(3 are integers are reported in the main body of the paper. 
Another relation of interest is the large \z\ expression 



$(a; /3; z) (-z)' 



r(/3- 



e z z a-0 



(A.6) 



To get this formula, we take the integral representation 
of $ 

$(a; /3; z) 

2m]_ ioo r(a) r(/3 + s ) i( Sjt Zj ' 

(A.7) 

where the poles of r(a + s) are on the left of the inte- 
gration contour and those of T(—s) are on its right. If 
\z\ < 1, the contour is closed to the right and one gets 
Eq. (|A.2j) . which is then analytically continued to all z. 
If \z\ > 1, one can close the contour to the left, and com- 
pute the integral by taking the residues at the poles of 
T(a + s): 

Res[r(s + a), s = -a - k] = ( , / , keN. 



k\ 



The result is 

$(a; /3; z) 



A I 



(-*)- E 



fe=0 



r(a + fc) 



r(/3) 



TfjS-a-fc) fc! 



, (A.8) 



which would not converge if M = +oo. In practice, one 
truncates the series up to some finite order M. The first 
contribution in Eq. (|A.6|) is the leading term k = 0. The 
second (subleading when z < 0, dominant otherwise) can 
be inferred via Rummer's transformation 



$(q; /3; z) = e z <&(/3 - a; 13; -z) . 



(A.9) 



giving 



M 



x (-^)' fc fA10 ) 

where we have truncated the series to order M. 

Combining these formulae in Eq. (|A.3p and using the 
identity 

T(l-a + k) , ,. k T(a) 



r(l-a) 



= (-ir 



r(a - k) 



one also gets 



*(a;/3;z) £>) fc (l + a - (3) k ( " ' 

k=Q 



k\ 



(All) 

The phases of Eqs. ([76|) and (|7T[) can be inferred as fol- 
lows. One chooses the phase of z near the branch cut 
of W and the position of the latter. Let them be e l7r 
and the negative real axis. From the above equation, 
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02 ~ t^ 2 z e / 4 Z-^ 4 = i«/3 e »*rfl/4 When z > 0, n = 0, 

while when z < the phases of z e / 4 and z~ 6 / 4 have to 
be taken above and below the branch cut, so that n = 2. 

Formulae 13.4.10 and 13.4.23 of [97[ are helpful to write 
the equations of motion in a simple form: 

(a + zd z )®(a; /?; z) = a$(a + 1; j3; z) , (A.12) 
(a + zd z )^{a; (3; z) = a(l + a - /3)*(a + 1; /?; z) . 

(A.13) 

The derivative of with respect to its first argument is 
used to get the solutions in a four-dimensional power-law 
universe: 

d a Ha; M = J^y X>(° + k) - tf(a)] j^-^ ^ , 

(A. 14) 



where "0 is the digamma function and 



^{v + k)-^{v)=Y j — -. (A.15) 



When u = l, H k = tp(l + k) — 0(1) = 1 + 1/2 + 1/3 + 
• ■ • + 1/fc is called the harmonic number. 

Another function used in this paper is the generalized 
hypergeometric function 

t? ( a x \ - ( a )k(P)k z k 

,fi(o,ftT,,i ' )a ij(W«' ( 5 
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